Water-wave gap solitons: An approximate theory and accurate numerical experiments 
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It is demonstrated that a standard coupled-mode theory can successfully describe weakly- 
nonlinear gravity water waves in Bragg resonance with a periodic one-dimensional topography. 
Analytical solutions for gap solitons provided by this theory are in a reasonable agreement with 
accurate numerical simulations of exact equations of motion for ideal planar potential free-surface 
flows, even for strongly nonlinear waves. In numerical experiments, self-localized groups of nearly 
standing water waves can exist up to hundreds of wave periods. Generalizations of the model to the 
three-dimensional case are also derived. 

PACS numbers: 47.15.K-, 47.35.Bb, 47.35.Lf 



I. INTRODUCTION 

As we know from the nonlinear optics, specific self- 
localized waves can propagate in periodic nonlinear me- 
dia, with a frequency inside a spectrum gap. These waves 
are referred to as gap solitons (alternatively called Brag g 
solitons; see, e.g., Refe.0, H |, ! 1, |, 0,1, I, M EJ). 
Also in the field of Bose-Einstein condensation, gap soli- 
tons (GS) have been known [l2], [H, [HI- Recently, it 
has been realized that GS are also possible in water- 
wave systems [l5j . In particular, very accurate numeri- 
cal experiments have shown that finite-amplitude stand- 
ing waves over a periodic one-dimensional topography 
are subjected to a modulational instability which spon- 
taneously produces Bragg quasisolitons — localized co- 
herent structures existing for dozens of wave periods. 
However, in the cited work [15j . no analytical approach 
was presented. As a result, many important questions 
about water-wave GS were not answered, concerning 
their shape and stability. The present work is intended 
to clarify this issue, at least partly. More specifically, for 
a given periodic bottom profile with a spatial period A, 
we shall derive, in some approximation, coefficients for a 
standard model system of two coupled equations, describ- 
ing evolution of the forward- and backward-propagating 
wave envelopes A±(x,t) (see, e. g., Refs.pl. \4. Ho|b 

i(d t ±V g d x )A± = AA T + (T s \A ± \ 2 + r x \A T \ 2 )A ± , (1) 

where t is the time, x is the horizontal coordinate in the 
flow plane, and A±(x,t) are slow functions. Let at equi- 
librium the free surface be at y = 0. Then elevation of 
the surface y = rj(x, t) is given by the following formula, 



r](x,t) = Re 4 +e «-^(^)' + 4_ e -'"^oW< 
+ higher-order terms in kA±, 



(2) 



where n = 2-7r/(2A) is the wave number corresponding to 
the main Bragg resonance, ll)q(k) = [c//t tanh^QK)] 1 / 2 is 



the frequency at the gap center, g is the gravity accel- 
eration, and ho is an effective depth of the water canal 
[definitely, ho is not a mean depth; more precisely it will 
be specified later by Eqs.([6]) and (J7J]. The coefficients in 
Eqs.JT]) are: an effective group velocity V g = oIluq(k) / dn, 
a half-width A of the frequency gap, a nonlinear self- 
interaction r<j, and a nonlinear cross- interaction ■ 

Generally, it is assumed in derivation of the above sim- 
plified standard model that: (a) dissipative processes are 
negligible, (b) a periodic inhomogeneity is relatively weak 
(that is A <C Wo), (c) the waves are weakly nonlinear, 
(d) original (without inhomogeneity) equations of mo- 
tion, when written in terms of normal complex variables 
At, contain nonlincaritics starting from the order three: 

iA k m w*(kL4 k + i J T(k,k 2 ;k 3 ,k i )Al 2 A k3 A ki 

xS(k + k 2 -k 3 -k4)dk 2 dk 3 dk 4 , (3) 

where w* (k) is a linear dispersion relation in the absence 
of periodic inhomogeneity (a weak inhomogeneity adds 
some small terms to the right hand side of Eq.([3]); the 
most important effect arises from a term LA k , where 
L is a "small" linear non-diagonal operator). It is also 
required, (e) the coefficient T(ki, k%; k±) of the four- 
wave nonlinear interaction should be a continuous func- 
tion. In application to water waves the requirements (d) 
and (e) mean that: (i) all the second-order nonlincar- 
ities are assumed to be excluded by a suitable canoni- 
cal transformation (the corresponding procedure is de- 
scribed, e.g., in Refs.[H, [13); (ii) the model system ([T]) 
can be good only in the limit of relatively deep water, 
since on a finite depth the function T(ki, k2; k3, kj) is 
known to contain discontinuities which disappear on the 
infinite depth (see, e.g., Ref.[I3|)- Therefore we introduce 
a small parameter 



e = exp(— 2nho) -C 1, 



(4) 
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and we consider in the main approximation only the 
principal effect of weak spatial periodicity, namely cre- 
ation of a narrow frequency gap with A ~ eloq under 
the main Bragg resonance conditions. Then we imply 
a standard procedure for obtaining approximate equa- 
tions for slow wave envelopes, where the deep-water limit 
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of T(kx, k^; IC3, ki) is used for the coefficients T$ oc 
T(k, k; k, k) and Tjf oc 2T(k, — n; k, — k). Thus we ne- 
glect in the actual nonlinear wave interaction some rela- 
tively small terms with coefficients of order e. 

Of course, the functions A± should be sufficiently "nar- 
row" in the Fourier space, since dispersive terms propor- 
tional to second-order derivatives d^.A± are not included 
into the model. 

After derivation of all the coefficients in section II, 
some known "solitonic" solutions of Eqs.([T]) will be com- 
pared to numerical results for exact hydrodynamic equa- 
tions, with nearly the same initial conditions as in the 
solitons (in section III) . We shall see that very long-lived 
self-localized groups of standing water waves are possible. 
In some region of soliton parameters, water-wave GS ex- 
ist up to hundreds of wave periods, until unaccounted 
by Egs.([T]) processes change them significantly. In sec- 
tion IV we discuss some promising directions of further 
research, concerning three-dimensional generalizations of 
the coupled mode equations. Some auxiliary calculations 
are placed in two Appendices. 



II. COEFFICIENTS OF THE MODEL 

We start our consideration with a short discussion 
of conditions when dissipation due to bottom friction, 
caused by water (kinematic) viscosity u, is not impor- 
tant in wave dynamics. Obviously, a viscous sub-layer 
should be relatively thin in this case: db -C A. In a 
nearly linear regime, a width of the sub-layer can be esti- 
mated as db ~ (z^/cj) 1 / 2 , where u ~ (g/A) 1 / 2 . This gives 
us the following necessary condition for applicability of 
the conservative theory: 



A 3/4 1/4^-1/2 ^ L 



(5) 



Generally speaking, one cannot exclude a possibility that 
in a strongly nonlinear regime the vorticity can some- 
times be advected by a wave-produced alternating veloc- 
ity field far away from the rigid bottom boundary. Such 
vortex structures are typically generated near curved 
parts of the bed, and they can significantly interact 
with surface waves. However, we assume this is not the 
case; otherwise, the problem becomes too complicated. 
Though we do not have simple criterion to evaluate influ- 
ence of the bottom-produced vorticity, with A > I m we 
still hope to be correct when neglecting water viscosity, as 
well as compressibility and surface tension. This allows 
us to exploit the model of purely potential free-surface 
ideal fluid flows, commonly used in the water wave the- 
ory. 

Since in this work we consider the case of relatively 
deep water, we can write u>q(k) sw w*(/t)(l — e), where 
w*(k) = (gft) 1 / 2 is the frequency corresponding to the 
infinite depth. Later we will see that values e = 
0.01 . . . 0.02 are of the most interest. 

Let us introduce conformal curvilinear coordinates 
(C11C2) determined by an analytic function B(Q, with 



C = Ci + i((2 - ha), so that 

00 

x + iy = B(C) = C - 2k" 1 J2 l 3 ^ 1 sin(2n K C), (6) 

n=l 

with real coefficients f3 n . Without loss of generality, we 
assume /3i > 0. The unperturbed water surface y = 
corresponds to real values of C = Ci ~~ *0, while at the 
bottom we have (2 = 0, and 



xW(Ci) + zyW(Ci)=B(C 1 -^o) 



(7) 



is a parametric representation of the bed profile, which 
can be highly undulating (see, for example, FigfTJ) . 

In these conformal coordinates, a spectrum of lin- 
ear potential waves is determined through the following 
equation (compare to Ref.[18j|, where an analogous ap- 
proach but slightly different notations were used): 



gkt&nh(h k) = 0, 



with 



B'(&) = 1 - 4]T n/3„e"cos(2n< 1 ). 



(8) 



(9) 



Here [k tanh(Ziofc)] is a linear operator which is diago- 
nal in Fourier representation: for any function /(Ci) = 
J fk exp(ik£i)dk/2Tr we have 



[fctanh(M)]/(Ci) = J ktaTLh(h Q k)f k e ikCl dk/27T. (10) 

The eigenfunction ^i(C0 takes the following form, 

+00 

*„(&)= J2 c « e2mKCl > ( n ) 



with some coefficients c„. With a given /x, we have 
an infinite homogeneous linear system of equations for 
c„. Non-trivial solutions exist for some discrete values 
<x>( m )(^). The first gap in the spectrum is the differ- 
ence between the two first eigenvalues at fi — k, that is 
2A = u!(2)(k) — W(x)(/t). Approximately, for small e these 
eigenvalues are determined by the coefficient j3± (compare 
to Ref.pJ]), 



w?j 2 )( K ) ~ <?Ktanh(/iQK)(l =F 2e(i\). 



(12) 



It should be noted that u>m(K) corresponds to ^ K 

sin(KCi): while W( 2 )( K ) corresponds to w cos(k£i). 
Thus, in the first order on e, the half-width A of the gap 
in the spectrum of linear waves is 

A w w*(k)e0i = w*(k)A, (13) 

where A = e/3i <C 1 is a small dimensionless quantity. 
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FIG. 1: Example I: unperturbed free surface and the bottom 
profile for 2h n = 1.4tt, D = 0.95, and eC = 0.01229 [see 
Eq.JH)]. 
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FIG. 2: Example I: free surface elevation at x = (at the 
center of GS) for 5 = 0.4. 



As to the nonlinearity coefficients Ts and T^, their 
values for the case of infinite depth (in other words, their 
zcroth-order approximations in e) can be easily extracted 
from Ref . [l9( : 



-w*(k)k 2 . (14) 
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FIG. 3: Example I: free surface profiles for different time mo- 
ments when elevation at x = is at maximum. 



It should be noted, the cited work [lj| relies on results ob- 
tained earlier by Krasitskii [l6| , who calculated kernels of 
so-called reduced integrodifferential equation for weakly 
nonlinear surface water waves (see also the paper by Za- 
kharov (l7j . and references therein). It is important in 
many aspects that for deep-water waves the coefficients 
Tg and Tx have the opposite signs, and their ratio is 
Ts/Tx « -1/2. 

With the same zeroth-order accuracy, the group veloc- 
ity is 



1 w*(k) 



(15) 



Now all the coefficients have been derived, and the sim- 
plified coupled-mode equations for relatively deep water 
waves in Bragg resonance with a periodic bottom take 
the following explicit form: 

, I ^ + ^) a + = Ac- + \ (|a+| 2 - 2|a_| 2 ) o+, (16) 



to* 
dt_ 

a;* 



2k J 

dx 
2 k 



Aa + + i(M 2 -2| a+ | 2 ) o_,(17) 



where a±(x,t) = KA±(x,t) are dimensionless wave am- 
plitudes. Analytical solutions are known for the above 
system (see @ II 0, US), describing moving localized 
structures, the gap solitons. In the simplest case the ve- 
locity of GS is zero, and the solutions essentially depend 
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FIG. 4: Example II: free surface elevation at x = for 2hoK 
1.2tt, D = 0.7, eC = 0.022, and 8 = 0.0. 



on a parameter 8, a relative frequency inside the gap 
(-!<£< 1): 



a± = y I(x) exp[— iSAuj^t + 170 ± i<p(x)], (18) 
4A(1 - S 2 ) 



I(x) 



<p{x) = arctan 



cosh[4A\/l - (5 2 



1 - 5 
T+6 



tanh 1 



(19) 
.(20) 



These expressions correspond to purely standing, spa- 
tially localized waves with frequency uj = (1 — e + 5A)lu* 
(concerning their stability, see Ref. 0, where, however, 
stability domains were presented for a different ratio 
Ts/rx; there are some numerical indications that the 
above GS are stable in a parametric interval <$* < 6 < 1, 
where a critical value (5* w —0.4). 

It should be noted, one can hardly expect a detailed 
correspondence between the very simple model (|16H17p 
and the fully nonlinear dynamics, but just a general ac- 
cordance sometimes is possible. In particular, the model 
does not describe nonlinear processes resulting in gener- 
ation of short waves which take the wave energy away 
from a soliton, thus influencing its dynamics. The model 
is also not generally good to study collisions between soli- 
tons, since wave amplitude can significantly increase in 
intermediate states. 
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FIG. 5: Example II: free surface profiles for different time 
moments corresponding to maximum elevation at x = 0. 
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FIG. 6: Example II: formation of sharp wave crests over bar- 
riers. 



III. NUMERICAL EXPERIMENTS 

In order co compare the above approximate analytical 
solutions to nearly exact numerical solutions, we chose 
the following function B(Q: 

with real parameters < Do < 1, and < C < 1. Hence, 
A = eD C < e. For D and C both close to 1, Eq.® 
gives periodically arranged barriers (see, for example, 
FigQJ. The barriers are relatively thin as (1 — Dq) <C 1, 
and relatively high as C — > 1. However, in numerical ex- 
periments with high-amplitude waves, a strong tendency 



was noticed towards formation of sharp wave crests over 
very thin barriers (say, when Do = 0.99), already af- 
ter a few wave periods. With sharp crests, the conser- 
vative potential-flow-based model fails (it is also clear 
that tops of narrow barriers must generate strong vortex 
structures). Therefore we took Do = (0-7 . . . 0.9) in most 
of our computations in order to have a smooth surface 
for a longer time. Exact equations for ideal potential 
free-surface planar flows were simulated (their derivation 
can be found in Ref . (l8| . so me g eneralizations are made 
in Refs.[2(| UM)- As in Ref. [HI], we dealt with dimcn- 
sionless variables corresponding to g = 1, K = 100. The 
dimensionless time t is then related to the physical time 
t = rt by a factor r = (lOOA/Vg) 1 / 2 . For instance, the 
period of linear deep-water waves with the length A = 2A 
is T* = (2tt/\/I00)t ~ 0.628r. At t = 0, we set the hor- 
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FIG. 7: Example III: interaction of two water-wave GS. 

izontal free surface, while the initial distribution of the 
surface-value velocity potential was 



MCi) = 2«T 3 / 2 /T(C 1 )cos[ K Ci + ¥>(£)] « 

(22) 

in accordance with approximate relation r\ t ~ «^>. 

Many simulations with different parameters were 
performed, and a very good general agreement was 
found between numerical and analytical results in the 
weakly-nonlinear case, that is for small steepness s = 
2[J(0)] 1 / 2 = 4[A(l-5)] 1 /2 < 0.35. So, with 2h n = 1.4tt, 
D = 0.95, eC = 0.01229, and S = 0.4 (example 
I), some noticeable deviations from the purely-standing- 
wave regime were observed only after t/r > 120 (see 




FIG. 8: Example III: highly nonlinear wave near x = 0. Here 
the free surface profiles are presented from t/r — 35.16 to 
t/r = 35.80 with the time interval At/r = 0.04. The wave 
profiles, except the last one, are given vertically shifted for 
convenience. Also the bed shape is shown. 



FigsOEJ). In a real- world experiment it could be several 
minutes with A~lm. 

What is interesting, even for larger s, up to s ~ 0.48, 
GS can exist for dozens of wave periods. A numeri- 
cal example for such a relatively high-amplitude water- 
wave GS is presented in Figsf4][5l where 2h$K = 1.2ir, 
D = 0.7, eC = 0.022, and 8 = 0.0 (example II). In this 
simulation, there were 45 oscillations before sharp crests 
formation (see Figs [3] and [5]). As to a further evolution 
of such GS, only in a real- world experiment it will be 
possible to get reliable knowledge about it, since various 
dissipative processes come into play. 

Concerning water-wave GS with negative 8, their be- 
havior for 8 > —0.4 was found stable, while for 8 < —0.5 
the dynamics was unstable, and partial disintegration of 
GS was observed after a few tens of wave periods (not 
shown). However, in some numerical experiments, the 
life time of GS was limited by the above mentioned pro- 
cess of sharp crest formation rather than by their own 
instability in frame of the model (|16fll7|) . at least with 
s 0.4 (not shown). 

Finally, we would like to present an example of inter- 
action of two GS (example III). The bed parameters are 
2h K = 1.2-k, D = 0.7, eC = 0.023. Both solitons ini- 
tially had 8 — 0.4 and they were separated by a distance 
66A. At t = we set the horizontal free surface and 



^o(Ci) w ^gs(Ci - 33A) + * GS (Ci + 33A). 



(23) 



This numerical experiment also describes interaction of 
a single GS with a vertical wall at x = 0. Surface pro- 
files for several time moments are shown in Figs|7][H] We 
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see that in this example the interaction between GS is 
attractive. They collide and produce a highly nonlinear 
and short wave group near x = 0. 



IV. 3D GENERALIZATIONS AND DISCUSSION 



the first pair corresponding to wave vectors ±pi = 
±(7r/A)(l, 1), and with the second pair corresponding to 
±P2 = ±(7r/A)(— 1, 1). It is important that the absolute 
values are equal to each other: |pi| = | p>2 1 = v2k = x. 
Again we will assume e = exp(— 2xh) < 1. It is conve- 
nient to use new horizontal coordinates: 



In this work, coefficients of the standard model ((T|) 
were derived for water-wave GS in the approximation of 
relatively deep water. The frequency gap in this case 
is small (of order e) despite strong bed undulations. It 
seems that a more general situation of intermediate depth 
cannot be described by this basic model, since an in- 
teraction of the main wave with a long-scale flow ("ze- 
roth harmonics" ) is then essential and should be included 
into equations. At the formal level, this corresponds 
to the mentioned discontinuities of the four-wave ma- 
trix element T(ki, k.2; ks, k.4) on a finite depth. Actually, 
in a finite-depth dynamics, three-wave interactions are 
more essential, and therefore they cannot be removed ef- 
ficiently by a weakly-nonlinear transformation. This is 
the main difference between the present third-order the- 
ory and previouslydeveloped second-order theories (see, 
for example, Ref. [221] ) . 

So far we considered purely two-dimensional flows, 
with the single horizontal coordinate x. Let us now intro- 
duce two important generalizations for three-dimensional 
flows. Below we only derive equations, but their detailed 
analysis will be a subject of future work. 

In the first case, the bottom topography is still one- 
dimensional, but we take into account weak variations of 
the wave field along the second horizontal coordinate q, 
simply by adding dispersive terms, proportional to <9 2 a±, 
to the coupled-mode system, as written below: 



id t , id x dg \ - 
— ± — + — % a± = A a 
w, 2k Ak 2 I 



i t - (|a± 



2 -2|a T | 2 ) a±. 



(24) 

In this system, a near-band-edge approximation for the 
upper branch of the linear spectrum gives a 2D focus- 
ing nonlinear Schrocdinger equation (NLSE). Thus, in a 
long-scale limit, the system (|24p exhibits a tendency to- 
wards wave collapse which is known as a typical feature 
of 2D NLSE dynamics. 

In the second periodic bottom profile 

Y^(x, q) = —h + \(x, q) is essentially two-dimensional, 
and in the horizontal Fourier-plane there are several 
pairs of Bragg-resonant wave vectors. For simplicity, we 
present below equations for the case when x(x, q) has 
the symmetry of a square lattice, with equal periods A 
in both horizontal directions x and q: 



otnin-z [cos(2niKx+2n2Kq)+cos(2niKx—2ri2Kq)] 1 

(25) 

where coefficients possess the symmetry a nin2 = a„ 2nr 
Let us consider interaction of two wave pairs having 
slow complex amplitudes a±(x,q,t) and b±(x,q,t), with 



V2 ' 



X2 



q — x 



(26) 



Elevation y = rj(xi, £2, t) of the free surface is then given 
by the formula 



xrj = Re^ e 



} + ..., 



6_e" 



(27) 



where Qq — \gn tanh(/i >*)] 1//2 ~ f2»(l — e ), with f7„ = 
(gx) 1 / 2 and eo = exp(— 2xho), and the dots correspond 
to higher-order terms (again, we should note that gener- 
ally ho 7^ h). Approximate equations of motion for the 
amplitudes have the following form: 



dt_ 


± 


d Xl 




2x 




± 


d X2 




2x 



a± = ei a T + ( 2 [b+ + bJ\ 



dHnl 

da* ± 
dHm 
dbX 



(28) 



(29) 



where small constants ei and £2 depend on a given 
bed profile, a|. and 6j_ mean the complex conjugate 
quantities, and the function H. n \ corresponds to non- 
linear interactions. Using an explicit expression from 
Ref. [13] for the deep-water four- wave resonant interac- 
tion T(ki, k2; k3, kj), we have 

n ni = I{|a + | 4 + M 4 + |6 + | 4 + |M 4 } 



-\a + \'\a^-\b, 

l 2 IM 2 - 



a , 



l 2 |6_l 2 



-T±{\ 

3 

-— [a+a-b*i_b*_ + a^alft+fe-] 



|a_| 2 |& + | 2 + |a_| 2 |6_| 2 } 



(30) 



where t± = Ti2i2(0) « 0.02346 is a normalized value 
of the matrix element T(ki, k2; ki, k2) for two perpen- 
dicular wave vectors of equal length (see Fig[3]). Since 
t_l <C 1, we actually may neglect in H. n i the terms pro- 
portional to T_|_. 

Unfortunately, it is hardly possible to find some ana- 
lytical space-dependent solutions for the nonlinear sys- 
tem (|28|) - (f29| . but it can be investigated by approximate 
methods. 

The parameters eoi £ii an d (-2 can in principle be cal- 
culated from solution of a linearized problem for water 
waves over a periodic 2D bed. An exact linearized equa- 
tion for a surface value of the velocity potential can be 
written in the following form: 



{u?jg - [fetanh(M)] -iV)^(r) = 0, 



(31) 
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FIG. 9: A normalized matrix element T(ki, k2; ki, k2) for two 
wave vectors of equal length, with an angle between them. 



where r = (x, q) is a radius- vector in the horizontal plane, 
k = (k%, + k^) 1 / 2 , while TV is a self-conjugate linear opera- 
tor corresponding to a bottom inhomogeneity. However, 
in three dimensions there is no compact form for N, valid 
with any bottom profile. At the moment, there are only 
approximate expressions N w Ni + N 2 + ■ ■ ■ + N m , ob- 
tained by expansion (up to a finite order m) of a vertical 
velocity at the level y = in powers of (relatively small) 
bottom deviation x( x ->q) from a constant level y = —h. 
The linear self-conjugate operators Nj have a general 
structure 



with 



Nj = [cosh(/iA:)]- 1 5' J [cosh(/ifc)]- 



Si = (V X V), 



S 2 = -(V X V) 



tanh(/ifc) 



Sa = (V X V) 



tanh(/ifc) 



(V%V) 



(VxV), 



tanh(Zifc) 



;(Vx 2 V)(V x V) - -(V X 3 V)V 
! 6 



(32) 

(33) 
(34) 

(V X V) 
(35) 



and so on, where V is the horizontal gradient (see Ap- 
pendix A). 

Now we are going to calculate £q, ei, and e 2 . Let us 
note that with e -C 1 the four independent eigenfunc- 
tions in Bragg resonance are: ^ cc w cos(kx) cos(nq), 
^ ss « sin(«;a;) sin(K(7), \l/ cs w cos(ra) sin(/cg), and ^ sc « 
sin(Kcc) cos(«;g). Accordingly, we have for the eigenfre- 
quencies 



ulJg sa ^tanh(/i>r) + (* cc 7V^ cc )/(*2 



(36) 



(where 



mean the average value in the horizontal 



plane), and analogously for w ss and u%. 



Let us 



introduce short notations for small quantities: v cc = 
(^ cc N^ cc )/[2x(^l c )} < 1, and similarly for v ss and 
fcs = ^sc- Then we have approximate equalities, 



n*(i 
n*(i 
n*(i 



(37) 
(38) 
(39) 



These frequencies should be identified with the eigenfre- 
quencies of the linear part of system (|28[) - (|29)h for space- 
independent solutions: 

1,1,1) = - eo) + ei + 2e 2 ] = w cc , (40) 

W(i,i,-i,-i) = 0*[(1 - e ) + ei - 2e 2 ] = w ss , (41) 

^(1,-1,1,-1) = - eo) - ei] = Wcsj ( 42 ) 

^(1,-1,-1,1) = - e ) - ei] = w sc . (43) 

As the result, we obtain the required formulas for the 
model parameters: 



1 



ei 



C2 



(^c 

('A!: 



(^cc + ^ss + 2f cs ), 



2za 



and 



(44) 
(45) 
(46) 
'cm is 



With Eqs.([211)-(IS3, calculation of v cc , 
straightforward if the function x( x il) contains a finite 
number of Fourier harmonics, for example 



X 



Oil 



«2 



cos 2kx + cos 2nq) H cos 2kx cos 2nq. (47) 



Moreover, since e <§; 1, it is possible to simplify the op- 
erators Sj by writing there instead of [tanh(/ifc)/fc]. 
By doing so and taking into account only Ni and N 2 , for 
the bottom profile (|47|) we obtain approximately 



(-2 



ei 



2e 

7! 



aia 2 , 



V2 



«2, 



2 v / 5 



(48) 
(49) 

(50) 



Thus, the expansion N « JVj + iVjj + . . . + TV m is cer- 
tainly useful for analysis of the case |V X | 1, but it can 
hardly be valid for a strongly undulating bed. It should 
be noted that a global representation of the velocity po- 
tential in the form ([54| (see Appendix A) is questionable 
in the general case. Derivation of N for arbitrary |V6|, 
assuming e -C 1, is an interesting open problem. 

It is worth noting that an explicit (though approxi- 
mate) form of operator N allows us to derive weakly 
nonlinear equations of motion for water waves over a 



8 



nonuniform 2D bottom. For example, the Hamiltonian 
functional (it is the kinetic energy JC plus the potential 
energy (g/2) J rfcPr) up to the 4th order in terms of the 
canonically conjugate variables ij(r 7 1) and tp(r, t) is writ- 
ten below: 



H 



ipKip + grj + 7] 



(W) 2 - (k^f] } 



d 2 r 



^Knknkijj + r] 2 (kip)V 2 ip d 2 r, (51) 



where K = [kt&nh(hk) + N] (see Appendix B). It is in- 
teresting to note that the bottom inhomogeneity comes 
into the Hamiltonian through the definition of opera- 
tor K only. For iV" = 0, it coincides with the previ- 
ously known fourth-order Hamiltonian for water waves 
on a uniform depth (see, e.g., Ref . [T^ . and references 
therein). It is also clear that coupled- mode system (|28|) - 
(|2"9")1 corresponds to the case K w f t, when the difference 
[K~ k] is neglected in the third- and fourth-order parts of 
the Hamiltonian, but it is kept in the second-order part. 
The functional Tl{ri,ip} determines canonical equations 
of motion, 



6H 



m = — ~ Ktj) 

+ Kr/Kr/Kip 



1 r 
2 



Kri 2 V 2 il> + V 2 rj 2 KiP 



,(52) 



5H 
8r) 



911 + 2 



(v^) 2 - (k^) 2 
+ j](k^)v 2 ij. 



(53) 



Numerical simulation of these cubically nonlinear equa- 
tions, with N 7^ 0, will be an important subject of future 
research. 

Further analytical and computational work is also 
needed to investigate formation of vortex structures near 
the bottom boundary and to evaluate their influence on 
the free surface dynamics. In any case, the present re- 
sults, based on the 2D purely potential theory, are deserv- 
ing much attention. Moreover, the author hopes that in 
a future real-world experiment all the mentioned dissi- 
pative processes will not be able to destroy water-wave 
GS for a sufficiently long time. Instead, with vortices 
and breaking wave crests, the predicted phenomenon of 
standing self-localized water waves over a periodic bed 
will be found even more rich, interesting, and beautiful. 
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Appendix A. Expansion of operator N 



The expansion of N in powers of x is easily obtained from the integral representation of the velocity potential 

d 2 k 



$(r,y) 



^ cosh A: (y + h) ^ sinhfcy 
coshkh k 



(2tt) 2 



(54) 



where 0k is the Fourier transform of the velocity potential at y = 0, and /k is the Fourier transform of an unknown 
function /(r) which should be determined by substitution of Eq. (|54[) into the bottom boundary condition 



ld<z>/dy-v x -v$}\ v= _ h+x 

The resulting integral equation can be represented as follows, 



= 0. 



, ikcosh[fc(x(r) - h)] ... . d 2 k 
/k P 6Xp(ikr W 



iksinh[fcx(r)l d 2 \t 
?k — rrm exp(zkr) 



k cosh[kh] 



0. 



(55) 



(56) 



It can be formally solved for /(r) by expanding Eq. (|56p in powers of x an d assuming / = fx + f% + For instance, 

equation (|56p with the third-order accuracy is written below: 



1 + (VxV) 



tanh(/ifc) 



- I V y v ) ) [coshfa*)]/ 



(VxV) - ( V^V ) V 2 



[cosh(ftA)]-V- 



(57) 



As the result, be obtain an approximate solution / w (Nx +N2 + ^3)0, where the operators Nx, N2, and N3 are given 
by Eqs. ((32)1 - (135)) . A linearized system describing the free-surface dynamics is 

- if>i = gr), rj t = [fctanh(M)]V> + /, (58) 

where ip(T) = $(r, rj(r)) is a surface value of the velocity potential (in the linear approximation tp w 0). It gives us 
the equation (f3"Tj) for eigenfunctions 9 U corresponding to some fixed frequency uj. 
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Appendix B. Hamiltonian of water waves up to 5th order 



An approximate Hamiltonian of water waves can be easily derived by writing the kinetic energy of potential three- 
dimensional motion of an ideal fluid in the following form, 

?)(r) 

fC = \ J d 2 r J [(d$/dy) 2 + (V<&) 2 ] dy = \ J ^ [d$/dy - Vr, ■ V<&] \ y=r] d 2 r 



= 5 /w- 



^ iksinh[fe(ry(r) + h)] ^ , ikcosh[/c77(r)] 
fccosh[fc/i] 1 k 2 

2 3 

ipk<j> + 77V tp ■ V(f> + — vip ■ vk<t> - — Vip ■ v(v 2 ^ 

2 6 



d 2 k 

6Xp(zkr) (2^)2 ^ 



d 2 r, 



with subsequent substitution 

4> rs V - vKip + r}kr]k^ + yV 2 ^ - r/k (r]kr]kip + yV 2 ^j - y V 2 
The approximate equality ([60|) follows from an expansion of Eq. ([54| : ip « 1 + 77A — ^-V 2 — ^-V 2 ^ 



i/Kib + —V 2 kib. 
6 



(59) 



(60) 



Thus, 



„2 



^ V72 „,\ »? V72, 



V> - 7?A> + T]Kr]Kip + -Jr-V^V - rjK-qKip + -jrV^tp - -t-V'riKip + -5-V Kip 



2 

rf 2 

f/> — 77^^? + riKr\Ki\) + — V ip 



(3 



+ — ■ - 77^) - — VV> ■ V(V 2 V>) + . . . }d 2 r. (61) 



After simplifying, we obtain H = | J" |^^"^) 



W (5) 



rj 3 - 
~6 



+ (777 + ?7 



(Vi/>) 2 - (A'V') 2 ] } d 2 r + + + 



where 



(kip)V 2 K^ - c/w,/Y/,/w,/\ r - ^-(V 2 i/>) 2 - if (Kr]K^)V 2 ^ - ,J -{Kil))V 2 {r]kiP) 



(62) 



d 2 r. (63) 



In the same manner, it is also possible to derive the Hamilt onian with a higher-order accuracy. 
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